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ABSTRACT 

We present the first ~ of a series - study of the evolution of galaxies in compact groups 
over the past 3 Gyr. This paper focuses on the evolution of the nuclear activity and 
how it has been affected by the dense environment of the groups. Our analysis is 
based on the largest multiwavelength compact group sample to-date, containing com¬ 
plete ultraviolet-to-infrared (UV-to-IR) photometry for 1,770 isolated groups (7,417 
galaxies). We classihed the nuclear activity of the galaxies based on optical emission 
line and mid-infrared diagnostic methods, as well as using spectral energy distribu¬ 
tion htting. We observe a 15% increase on the number of the AGN-hosting late-type 
galaxies found in dynamically old groups, over the past 3 Gyr, accompanied by the 
corresponding decrease of their circumnuclear star formation. Gomparing our compact 
group results with those of local isolated field and interacting pair galaxies, we find no 
differences in the AGN at the same redshift range. Based on both optical and mid-lR 
colour classifications, we report the absence of Seyfert 1 nuclei and we attribute this 
to the low accretion rates, caused by the depletion of gas. We propose that the ob¬ 
served increase of LINER and Seyfert 2 nuclei (at low-z’s), in the early-type galaxies 
of the dynamically young groups, is due to the morphological transformation of lentic¬ 
ular into elliptical galaxies. Finally, we show that at any given stellar mass, galaxies 
found in dynamically old groups are more likely to host an AGN. Our findings suggest 
that the depletion of gas, due to past star formation and tidal stripping, is the major 
mechanism driving the evolution of the nuclear activity in compact groups of galaxies. 

Key words: Galaxies: nuclei - Galaxies: interactions - Galaxies: groups: general - 
Galaxies: Seyfert 


1 INTRODUCTION 


Active galactic nuclei (AGN) are considered some of the 
most dramatic and enigmatic phenomena in the evolution 
of galaxies. It is now believed that most galaxies in the Uni¬ 
verse ma y host a supper-m assive black hole (SMBH) in their 
nucleus ( Kormendvll2004l ~) . thus stressing the importance of 
understanding the key mechanisms creating, fuelling and 
quenching them. Despite the diversity of the observed nu¬ 
clear activities, there is currently a general consensus that 
they are all different expressions of the s ame phenomena . 
According to the AGN unification scheme jAntonuccill 19931 : 
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lUrrv fc Padovanlll995l l. this diversity depends on our view¬ 
ing angle. The SMBH is surrounded by an accretion disk, 
which in turn is also enclosed in an optically thick dusty 
torus. Depending on our position, we can either observe di¬ 
rectly towards the central broad line region (BLR), where 
we may observe a Type 1 source (a.k.a Seyfert 1; Syl), or 
the torus shall block our view towards the center, allowing 
us to observe only its reflection on the narrow line region 
(NLR) clouds, and therefore a Type 2 AGN (a.k.a Seyfert 
2; Sy2). As a consequence, the expected spectrum of a Sy2 is 
dominated by narrow lines (<1000 km s“^) of ionized met¬ 
als (such as [OIII]A5007 and [NIIJA6584), whereas Syls have 
both narrow and broad lines (^1000 km s“^). 

More recently, however, a number of studies suggested 
that Syl and Sy2 hosting galaxies display important in- 
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trinsic , as well as environmental differences. ISteffen et al.l 
||2003D showed that broad-line AGN are dominating the 
high X-ray luminosities, whereas narrow-line AGN the lower 
ones, (the so-called “Steffen effect”). To interpret these re¬ 
sults one should assume either a modification to the uni¬ 
fied scheme, to include X-ray luminosities, or the com¬ 
plete segregation of Syl and Sy2 sources as two differ¬ 
ent categories of nuclear activity. Some optical and in¬ 
frared studies revealed that Sy2 galaxies are more likely 
to display circumnuclear s tar formation than Syls (i.e. 


Dultzin-Hacvan et al.l Il994l : iMaiolino et al.l Il997l : IGu et al] 
2001^ . with more than 50% of them exhibiting nuclear star- 


bursts dCid Fernandes et al.l [iooH i. Galaxies hosting a Sy2 
nucleus, though, were found more frequently in interac¬ 
tion with close neighbours (^100 kpc) than Syls or non¬ 
activ e galaxies i Dultzin-Hacvan et al. ]j99|;^Kron gojd et al.l 


I 2 OO 2 I : iKoulouridiset^jh 2006 : V illarroel fc KorrJ 2014 1. On 


the other hand. IWu et al.l ( 2009l i showed that there is 
statistical difference in the strength of the polycyclic aro¬ 
matic hydrocarbons (PAHs), which pr obe sta r form ation, 
between Syl’s and Sy2’s. In addition, lElitzuil l|2012ll pro¬ 
posed that Syl and Sy2 AGN are simply the extremes of 
the AGN distribution, where in the first case we observe 
directly towards the BLR where as in the oth er dust blocks 
completely our view towards it. iTranI ll200lll showed that 
50% of the Sy2’s in their sample displayed BLRs observed in 
reflected polarised light (a.k.a hidden-BLRs; HBLRs), con¬ 
firming their obscured Syl nature. However, the remain¬ 
ing non-HBLR galaxies appear to have lower A GN lumi¬ 
nosit ies, indicative of their lower accretion rates dWu et al.l 
[20 i 3). All the above results, raised concerns about the va- 
lidity of the Unification scheme. Alternative models (i.e. 
lNicastroll2000l : iNicastro et al.ll200.ll . and references therein) 
suggested that the accretion rate, and subsequently the 
AGN luminos i ty, pl ay a key role in the presence of the BLR. 
lElitzur et al.l (l200fJ i showed that for low bolometric lumi¬ 
nosities (<10'*^erg sec“^), the torus and possibly the BLR 
disappear, due to the inability of the SMBH to sustain the 
required cloud outflow rate. 


The existence of a connection between galaxy inter¬ 
actions and nuclear activity, both in terms of star for- 
mation and/or AGN ac t ivity, has been also confirmed by 
IStorchi-Bergmann et al.l (|2008h . Since most galaxies in the 
Universe are found in interac ting environment s - such as 
pairs, groups and clusters fe.g. ISmall et al.|ll999h - it is very 
important to understand their connection with the various 
types of nuclear activity. Although clusters host the largest 
number of galaxies, of all the interacting environments, it 
has been recently shown that a large number of their mem¬ 
bers have been pre-pro c essed in group s (i.e Gortese et al.l 
l2006l : lEckert et Sl2014l L lHickson et al.l (|l992 \ showed that 
compact groups (GGs) of galaxies appear as the ideal sys¬ 
tems to study the effects of galaxy interactions, due to their 
high galaxy densities and low velocity dispersions (-^250 
km s“^). So far, the most compl ete and best st udied CG 
sample was the one compiled by iHicksorJ (Il982l f . The ini¬ 
tial sample consisted of 100 GGs - the so-called Hickson 
com pact groups (HCGs), containing 451 galaxies - however 
later iHickson et ^ (ll992l L using spectroscopic information, 
reduced it to 92 groups having at least 3 accordant mem¬ 
bers. Although, the dynamical and star formation proper¬ 
ties of these groups have been extensively examined dur¬ 


ing the last three decades (i.e MgmjesjleOIiyeirafc HicksonI 
1 1994 Ijonhson et ahllioOTl : iBitsakis et al. I2OI0I . 2011^ . there 


is a li mited number o f stud ies focused on their nuclear ac¬ 
tivity. IShimada et al.l (l2000l ~) studied a sample of 69 galax¬ 
ies belonging to 31 HGGs and f o und t hat nearly 40% host 
an AGN. Later, iMartinez et al.l (120081 1 and iMartinez et al.l 
(l2010ll relied on the optical spectroscopy of 270 galax¬ 
ies in 64 GGs, and confirmed the fraction of AGN previ¬ 
ously found (42%). They also showed that the majority of 
AGN in GGs have low luminosities (LLAGN), and they at¬ 
tributed this to their hi gh gas deficiencies. This result i s 


consistent with those o f VerdeS::MQntenegro_et^ (2001), 
iMartinez-Badenes et al.l (l2012ll and Bitsakis et al. (2014), 
where it was shown that HGG galaxies display high deficien¬ 
cies in their atomic and molecular gas, as well a s their dust 
conten t due to tidal stripping. More recently, ISohn et al.l 
ll2013ll studied a sample of 58 local GGs selected from the 
SDSS data release 7. They found AGN activity in 17-42% of 
the galaxies, depending on the classification method used. 
They also showed that no powerful dust-obscured AGN were 


detected, accordin g to the most r ecent WISE mid-IR classi¬ 
fication method of iMateos et al.l (l2012li . and they suggested 
that nu clear activity is n ot st rong due to gas depletion. 
Finally, iTzanavaris et all (|2014h examined Chandra X-ray 
maps of 9 HGGs and showed that the majority of the galax¬ 
ies displaying low specific star formation rates (thus being 
earlier types), are also more likely to host a weak AGN, with 
Lx,o.5-8.ofceV < erg sec“^ 


Although these studies examined the incidence of the 
AGN activity in local GG galaxies (z<0.05) and its connec¬ 
tion with the dense environment of the groups, they were not 
able to show how it evolved throughout cosmic time. Never¬ 
theless, the well known con nection between star f ormation 
and the AGN activities fe.g. iMagorrian et al.lll9^ 'l and the 
fact that the star formation history of the Universe change d 
significantly over the past few Gyr (e.g. I Heavens et ahllioM L 
imply that AGN activity should have also experienced sev¬ 
eral different phases, during this period. Using the advent 
of wide-area extragalactic surveys, such as the Galaxy Evo¬ 
lution Explorer surveys (GALEX), the Sloan Digital Sky 
Survey {SDSS), the 2 micron all-sky survey {2MASS), 
and the Wide-field Infrared Survey {WISE) Explorer All 
Sky Release, as well as the published GG catalogues of 
iMcConnachie et al.l (l2009ll , we obtained a working sample of 
1,770 GGs, containing 7,417 galaxies, in the redshift range 
of 0.01<z<0.23 (thus, looking back time of approximately 
3 Gyr). Using this sample, we will examine the evolution of 
the AGN activity and its connection with the dynamic en¬ 
vironment of the groups. In §2, we describe the selection of 
our sample, and also the comparison samples we have used 
to put the properties of our galaxies into context. In §3, the 
nuclear and morphological classifications are presented, as 
well as the spectral energy distribution modeling we per¬ 
formed to estimate some important physical properties of 
the galaxies, necessary for our study. In §4, we present re¬ 
sults, and in §5 we discuss their implications. Finally, in §6, 
we sum up our main conclusions. 


To calculate the distances in this paper we adopt a flat 
ACDM cosmological model, with parameters: Ho=70 km s“^ 
Mpc“^, Dm =0.30 and Ha=0.70. 
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2 THE SAMPLES 

2.1 The SDSS compact group sample selection 

The sample of compact groups of galaxi es, presented in this 
paper is part of the larger sample of iMcConnachie et alJ 
(l2009l) . The latter was constructed by applying the slightly- 
improved Hickson’s criteria, to the wh ole Sloan Digital Sky 
Survey Data Release 6 (SDSS'-DR6; lAdelman-McCarth^ 
l2008h . iHicksonI lll982li defined as compact all galaxy groups 
in which, (i) the number of galaxies within 3 magnitudes 
are at least four, N(Amr=3)^4, (ii) the angular diameter 
with no additional galaxies is at least 3 times the diam¬ 
eter of the group, 6n -5 36g, and (Hi) the total magni¬ 
tude of the galaxie s, averaged over Og, is les s than 26mag 
arcsec”^. However. IMcConnachie et alJ (l2008l f showed that 
by increasing the required minimum surface brightness, from 
26-to-22 mags arcsec”^, the contamination of groups con¬ 
taining an additional non-member galaxy decreases from 
71% to 26%. Moreover, they imposed a faint-end limit at 
Petrosian r-band magnitude of 18 mags, to ensure that the 
observational and the mock catalogues (created to test the 
selection criteria) are robustly compared, and a bright-end 
limit of 14.5 mags, to ensure that automatic de-blending of 
galaxies is reliable. Their final catalogue comprised 2,297 
co mpact groups, containing 9,713 galaxies (see Catalogue A 
of IMcConnachie et al.l 120091 '). For the purposes of this work 
we have constructe d a su b-sample of the Catalogue A of 
IMcConnachie et ^ (l2009l i. using as a criterion the avail¬ 
ability of ultraviolet (UV), optical, near and mid-infrared 
(IR) photometry, which are essential to estimate the physi¬ 
cal properties of the galaxies and perform a detailed analysis 
of the sample. 

Initially, since Catalogue A originates from the SDSS, 
optical phot ometry was available for all the galaxies. The 
SDSS-T)R7 (lAbazaiian et al.]|2009l l covered over 35% of the 
sky in the u,g,r,i,z bands, centered at 3557, 4825, 6261, 
7672 and 9097%, with a resolution of 1.4 arcsec and sen¬ 
sitivities ranging from 22.0-to-20.5 mags for 5 -ct detection, 
respectively. The photometry was automatically estimated 
using a combination of two models, a pure de Vaucouleurs 
profile and a pure exponential profile, which is considered 
the opti mal for extended so urces (called the CModel Mag¬ 
nitudes; (Strauss^L^HooS)- Simultaneously, spectroscopy 
was performed at all sources brighter than 17.77 mags in r- 
band, with angular separation more than 55 arcsec to avoid 
collision of the SDSS fibers. The size of the SDSS aperture 
was 3” (see more in §3.4). The fluxes and equivalent widths 
of the optical lines have been obtained from the work of 
iBrinchmann et ^ ll2004h . These authors have published the 
processed line fluxes of more than 8T0® galaxies from the 
SDSS-DR7 (see MPA-JHU DRiQ). They fitted single Gaus¬ 
sian profiles after removing the stellar continuum, correct¬ 
ing for the foreground galactic reddening and re-normalising 
them to match the photometric fiber magnitudes in the r- 
band. In addition, they re-estimated the uncertainties to 
include both errors in spectrophotometry and continuum 
subtraction uncertainties. Unfortunately, they have chosen 
500 km s“^ as an upper limit of the width of the Gaussian 
profile, to ensure the correct estimation of the faintest line 

^ http://www.mpa-garching.mpg.de/SDSS/DR7 


fluxes where widths are often overestimated, thus prevent¬ 
ing us from classifying real broad-line galaxies (with widths 
>1000 km s“^). The fraction of the emission line galaxies in 
our sample, having signal-to-noise ratios S/N>3, was 75%. 

For the near-IR observations we have used the pub¬ 
lished magnitudes from the 2MASS extended source cat¬ 
alogue llSkrutskie et al.ll2006li . The galaxies were imaged in 
J, H, and K bands, with 3-a sensitivity limits of 14.7, 13.9, 
and 13.1 mags, respectively. The photometry was performed 
using sets of circular and ellipsoidal apertures applied on the 
3-cr J-band isophote of each individual galaxy, so that both 
the integrated flux, as well as the estimated uncertainty were 
estimated. _ 

Using the latest WISE ilWright et ahllioi^ all sky re¬ 
lease, we obtained the photometry or the upper limits for the 
galaxies in our sample. WISE was a 40cm telescope, which 
performed all sky survey, covering the 3.4, 4.6, 12.0, and 
22.0^m bands (named wl-to-w4), with 5-a sensitivity lim¬ 
its 19.70-to-14.40 mags, respectively. The WISE photom¬ 
etry (WPHOT) was automatically performed using profile 
fitting. However, if the PSFs of two or more galaxies were 
closer than the resolution of the telescope (6-to-12 arcsec), 
WPHOT reports the total photometry with a blending pa¬ 
rameter, which is the number of PSFs included in the aper¬ 
ture. 

The UV data were obtained from GALEX 
llMorrissev et al.l l2005li All Sky Survey (AIS; with lOOsec 
on-source exposures). Medium Imaging Survey (MIS; 
with 1500sec on-source exposures), Nearby Galaxy Survey 
(NGS; with typical exposures of 4000sec), Deep Imaging 
Survey (DIS; with typical exposures of 4000sec), as well as 
from Guest Investigator’s data l-to-4 (GIl-4; with typical 
exposures of 1500sec); publicly available in the GALEX 
archive. GALEX is a 50cm diameter UV telescope that 
imaged the sky simultaneously in both FUV and NUV 
channels, centered at 1540A and 2300A, respectively. 
The field-of-view (FOV) is approximately circular with a 
diameter of 1.2° and a resolution of about 5.5” (FWHM) 
in the NUV. The data sets used in this paper are based-on 
the GALEX sixth data release (GR6). The photometry 
was initially acquired from the GALEX archive, where 
it was automatically performed using SExtractor (see 
iBertin fc Arnout^l 19961 '). a code that estimates the fluxes of 
all the sources in a given image using ellipsoid apertures. 

To check the accuracy of the archival photometry ob¬ 
tained in the various bands (from the UV-to-IR), we ran¬ 
domly selected a sample of galaxies and performed aperture 
photometry, using the galaxy 3 -ct isophotal contours. Our 
results suggested that the published SDSS, 2MASS and 
WISE archival fluxes were consistent with our photometry, 
within the uncertainties. However, for more than 10% of the 
GALEX/N\JX fluxes and for ~34% of the GALEX/FEY 
fluxes we found differences of 1-2 magnitudes (see Fig. [T]). 
To overcome this limitation in the UV fluxes, we manually 
performed aperture photometry by calculating the isophotal 
contours around each source in order to account for varia¬ 
tions in the shape of the emitting region. Then we defined 
a limiting isophote 3 -ct above the local overall background, 
for each galaxy, and we measured the flux within this re¬ 
gion, after subtracting the corresponding sky. To convert 
from counts to UV fluxes we used the conversion coefficients, 
given in the header of each file. The upper limits depend on 
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Figure 1. Differences between the GALEX photometry of the 
galaxies in our sample, derived manually using aperture pho¬ 
tometry (indicated as phot) and the archival ones (indicated as 
GALEX) in the FUV (top panel) and NUV (bottom panel) 
bands. In the bottom boxes are presented the standard devia¬ 
tions (cr) of the values created as running averages, in bins of 0.1 
mag. (A coloured version of this figure is available in the online 
journal) 
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Figure 2. The BPT diagnostic diagram of 4,208 galaxies in our 
sample with available optical spectroscopy. T he da shed and dot¬ 
ted lin es indicate the iKauffmann et al.l i2003l l and iKewlev et ^ 
ll2006ll AGN selection criteria, respectively. Galaxies located un¬ 
der the dashed line are classified as star forming nuclei (SFN), 
galaxies over the dotted line as pure AGN, and galaxies between 
the two as composite objects. Contours are at 5% intervals of the 
maximum of the total distribution. Mean uncertainties are indi¬ 
cated in the bottom right corner and correspond to the 6% and 
16% of the X and y-axis values, respectively. (A coloured version 
of this figure is available in the online journal) 


the survey ranging from 20.5 and 21.5 mags, to 24.0 and 
24.5 mags, for FUV and NUV bands respectively. 

Our final sample comprises 1,770 groups, containing 
7,417 galaxies with available UV-to-mid-IR flux densities. 
Due to SDSS fiber collision constraints, optical spec¬ 
troscopy was available only for 4,208 of these galaxies (for 
the purposes of this paper, galaxies with no spectral infor¬ 
mation were assigned with the z of the most massive galaxy 
in their group). This sample is nearly 55 times larger than 
the previous HCG multiw avelength samples presented in 
iBitsakis et al.l (1201 ll . |2014| ) . It also covers a much larger vol¬ 
ume in the Universe, with groups being detected up to a 
redshift of 0.23, equal to a distance of ~1.15 Gpc (the cor¬ 
responding for HCGs is z~0.022 or 95 Mpc). 


2.2 The comparison samples 

As comparison samples we have used truly Isolated galax¬ 
ies in the field, as well as galaxies in isolated interacting 
pairs. The first sample contains 513 local (with z below 0.05) 
isolate d field galaxies presented in iHernandez-Ibarra et al.l 
(l2013t) . It mostly consists of late-type galaxies (86%), with 
an AGN fraction of 40%, at least for the 87% of them for 
which optical emission lines could be measured. The second, 
is a sample of 385 isolated interacting pair galaxies, 71% 
of which are LTGs (originally selec ted from the catalog of 
lKnr^hentsevl[l97^ . was taken from iHernandez-Ibarra et al.l 
1 20141) . Approximately 82% of these galaxies have emis¬ 
sion lines, and ~49% of them are classified as AGN-hosting 
galaxies. The authors separated the pairs into three subsam¬ 
ples depending on the morphologies of the member galaxies. 
Pairs consisting of late-type galaxies (S-bS), of early-type 
galaxies (E-l-E), and mixed pairs (S-l-E). Both samples were 


selected from the SDSS-DK7 survey, having the same spec¬ 
troscopic limitations than our compact group sample. 


3 DATA ANALYSIS AND CLASSIFICATIONS 
3.1 Morphological classification 

A very important step to put the properties of our galaxies 
into context, was the classification of their morphologies, as 
we ll as of the dynamic al state of the groups (as proposed 
by iBitsakls et al.l 120101 ) . To classify our galaxies into late- 
type (spirals and irregulars; hereafter LTGs) and early-type 
(lenticulars and ellipticals; hereafter ETGs), we hav e relied 
on the criteria described in ISimard et al.l (l2009l . [201ll) . These 
authors have fitted the SDSS r-band 2D, point-spread- 
function-convolved, bulge-bdisk decompositions of a sample 
of 1.12 million galaxies from the SDSS-DR7. Using four dif¬ 
ferent decomposition procedures they ensured more robust 
structural parameters of the galaxies, even in crowded envi¬ 
ronments. According to their findings, galaxies with bulge- 
to-disk ratios B/T^O.35 and image smoothness at half-light 
radius 82^0.075 are classified as ETGs. Applying this clas¬ 
sification in our sample, 3,045 galaxies are classified as LTGs 
(41%), 4,367 as ETGs (59%). There is no available classifica¬ 
tion for 5. The results are similar to what we have observed 
in HCGs (54% LTGs and 46 % ETGs, respe c tively ). 

As mentioned above, in IBitsakis et al.l ll201(]|) we have 
proposed a dynamical-evolution classification for the groups, 
depending on their ETG fraction. There we showed that this 
kind of classification is physical and also consistent with pre¬ 
vious classifications based o n the distribution of the HI gas 
content of the groups (see IVerdes-Montenegro et ahl 1200 ll : 
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Figure 3. The iKewlev et al.l <l2006ll diagnostic BPT diagrams, 
that use [Ol]A6300 (top panel) and [Sll](A6717 +6731) lines (bot¬ 
tom panel), to further separate AGN galaxies into Sy2 and LIN- 
ERs (found in the upper left and lower right portion of both 
diagrams, respectively). Contours are at 5% intervals of the max¬ 
imum of the total distribution. Mean uncertainties are indicated 
in the bottom right corner of each panel and correspond to the 
18% and 16%, as well as the 6% and 16% of the x and y-axis 
values of the top and bottom panels, respectively. (A coloured 
version of this figure is available in the online journal) 



iBorthakur et aI.]l201Cll ). According to this, groups with more 
than 25% of early-type members are classified as “dynam¬ 
ically old” (DO), whereas groups with less or no ETGs as 
“dynamically young” (DY). Applying this classification to 
our current sample, we find 373 dynamically young (21%), 
and 1,397 dynamically old groups (79%). 


3.2 AGN classification methods 


To characterise the nuclear activity of the galax¬ 
ies in our sample, we have relied on the spectral 
classification diagnostic diag r ams, initially proposed by 
[Baldwin . Phillips fc Terlevichl (I l98ll ) and later expanded 
by Veilleux fc Osterbrock ( 19871 ) . also known as BPT- 
diagrams. These, compare the flux ratios of several nar¬ 
row emission lines, such as [Oiii]A5007/H/3, [Nii]A6584/Ha, 
[Sll](A6717 -E6731)/Ha and [Ol]A6300/Ha, and distinguish 


Figure 4. The lMateos et ahl ||2012|) WISE [4.6]-[12] versus [3.4]- 
[4.6] colour diagram. Galaxies containing dust obscured AGN are 
located in the area within the dashed lines. Gontours are at 5% 
intervals of the maximum of the total distribution. Since, <1% of 
our sources are found in the AGN locus, we identify them with 
black points. Mean uncertainties are indicated in the bottom right 
corner and correspond to the 4% of the x and y-axis values. (A 
coloured version of this figure is available in the online journal) 


the nuclear activity of galaxies between star-forming nu¬ 
clei (SFN), active galactic nuclei (AGN), and compos¬ 
ite objects (objects showing both signatures of both SFN 
and AGN). In Fig. [5] we present the [Oiii]A5007/H/3 ver¬ 
sus the [Nll]A6584/Ha line ratios of the 4,208 galaxies 
in our sample w i th av ailable optical spectroscopy (see 
iBrinchmann et ^ |2004| ). Follow i ng th e sele ction criteria 
introd uced by Kauffmann et all (l2003l) and iKewlev et al.1 
ll2006l) . 1,436 of them are classified as SFN (34%), 1,860 as 
AGN (44%), and 912 as composite objects (22%; detected 
to host both star formation as well as AGN activity within 
our aperture). 

However, the above methods are not able to distin¬ 
guish between broad-line Syl, which appear to have very 
wide spectral lines with widths >1000 km s“^, form narrow- 
line Sy2 sources. Never theless, as we explained in §2.1, 
IBrinchmann et H] (l2004h set an upper limit (of 500 km s“^) 
to the width of the Gaussian profiles they have used to fit 
the spectral lines, to better constrain the weak narrow lines. 
Therefore, using these data, we have classified as “broad¬ 
line AGN”, all AGN-hosting galaxies detected with signal- 
to-noise ratios >10a, and lines wider than the upper limit, 
a >500 km s“^. This had as a result to include both Syls, as 
well as other narrower sources (having spectral line widths 
between 500 and 1000km s“^). Yet, even applying this ap¬ 
proximate classification, the fraction of “broad-line” AGN 
in our sample is negligible (<1%). 

More recent classifications were also able to distinguish 
between Sy2s and low-ionisation nuclear emission-line re¬ 
gions (LINERs; their class ification as AGN is still 
open to discussion; e.g. see lHeckman|[l980l : IHo et al.1 
Il993l) . In the two panel s of Fig. [3l we presen t the diagnos¬ 
tic methods proposed by IKewlev et al.1 (l2006l ). According to 
these, we plot the [Oiii]A5007/H/? versus the [Oi]A6300/Ha 
(top panel) as well the [Sii](A6717 + 6731)/Ha (bottom 
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Figure 5. Example of fits obtained with the CIGALE code. The 
color code of the data-point refers to the instrument used, in 
purple: UV data from GALEX, in blue: optical data from SDSS, 
in red: NIR data from 2MASS, and in green: MIR data from 
WISE. The black solid line shows the best fit model for each 
source. The value of the reduced indicated for each fit. (A 

coloured version of this figure is available in the online journal) 


panel) line ratios, of all the galaxies detected in the clas¬ 
sic BPT as AGN-hosting. Based on these classifications, 
319/531 of these galaxies are classified to host a Sy2 nu¬ 
cleus (8-13%), whereas 685/851 are LINERs (16-20%), re¬ 
spectively. 

Finally, by app lying the colour d i agnos tic methods of 
iMateos et al.l (|20ld) and IStern et al.l (|2012h we were also 
able to identify dust obscured AGN activity. These meth¬ 
ods compare the mid-IR colours of galaxies, to identify a 
power-law-based selection of luminous AGN candidates. The 
first method compares the WISE [4.6/rm]-[12/rm] versus the 
[3.4/rm]-[4.6/rm] colours (see Fig. |4| and defines an area 
where power-law dominant AGN should be located. The sec¬ 
ond, describes different colour selection, where mid-IR lumi¬ 
nous AGN should be located, at [3.4/im]-[4.6^m]^0.8 (Vega 
mags). Both methods predict the fraction of dust obscured 
AGN to be <1%. 


3.3 Spectral energy distribution modeling 

We performed a UV to mid-IR broad band spectral energy 
distribution (SED; see Fig. [S| fitting of all of the sources 
of the catalogue using the recently updated version of the 
CIGALE model, which includes a detail ed treatment of 
the AGN influence to the observed SED dNoll et al.ll2009l : 
ICiesla et al.ll2015l . Boquien et al. in prep, Burgarella et al. 
in prep). Based on an energy balance between the energy 
absorbed in UV-optical and re-emitted in IR, CIGALE 
builds SED models that are then compared to the data. 
For each model, CIGALE computes the value in order to 
build the probability distribution function (PDF) of each 
parameter. The final value of a parameter is thus the mean 
value of the PDF and the uncertainty associated is the 
standard deviation of the distribution. 

We used a delayed star formation history (SFH) for our 
galaxies as it has been shown that it provides a good esti¬ 
mate of the stellar mass and SFR, as well as realistic ages 
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Figure 6. Selection biases in our sample. Malmquist bias pre¬ 
vents us from observing less massive galaxies at higher red- 
shifts. Moreover, due to an upper limit cut-off in the luminosi¬ 
ties, during the selection process, brighter galaxies at lower-z’s 
were also excluded. Galaxies selected in the stellar mass range of 
10.4:Clog(Mstar)^11.3 Mq, are located within the dashed lines. 
For the purposes of comparison we separate our sample into four 
different redshift regions, containing a similar number of galaxies 
(~1,100). These are: Binl (blue; 0.010<z<0.078), Bin2 (green; 
0.078^z<0.104), Bin3 (orange; 0.104s;z<0.133) and Bin4 (red; 
0.133<z<0.230). (A coloured version of this figure is available in 
the online journal) 


ijCiesla et al.ll201Rh . This SFH was convolved with the stellar 
population models of Maraston _|200^}Vrhe SED was then 
attenuated using the lCalzetti et alT i ioOfll l law, and the en¬ 
ergy a bsorbed is re-emitted in the IR through the lPale et al.l 
(l2014h templates. CIG ALE models the pr esence of an AGN 
through the library of iFritz et al.l (|200fil 'l. which takes into 
account the emission of the central object, the dust torus 
emission, and the scattered photons. However, when using 
broad band photometry, the AGN models are h ighly de¬ 
generated. Therefore, following ICiesla et al.l (l2015l l. we used 
three models: a Syl, a Sy2, and an intermediate type AGN 
template, displaying a power law emission in mid-IR without 
any strong UV contribution. Two of the most important out¬ 
put parameters, for our current study, provided by GIG ALE 
are the stellar masses as well as the contribution of the AGN 
to the total IR luminosity, Jagn, for each galaxy. We are 
confident in the presence of an AGN when / agn is higher 
than 10%. Between 5 and 10% ICiesla et al.l (l2015l l showed 
that this confidence is reduced (yet the x^ value is still 50% 
lower than by not using the power-law template), and for 
fractions below 5% we cannot conclude by the fit whether 
an AGN is present. 


3.4 Bias control and mass selection 

Due to its large size, with galaxies found in a wide range 
of redshifts and stellar masses, our sample can be subject 
to selection biases. In an effort to detect and control them, 
we plot in Fig. [6] the stellar masses of the galaxies in our 
sample as a function of their redshift. The Malmquist effect 
appears at higher redshifts (see lower-portion of this figure). 
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Since, we are unable to detect the faintest sources and, 
as a consequence we are biased towards the more massive 
ones. Furthermore, on the upper left portion of the same fig¬ 
ure we can notice one more bias. As we mentioned in §2.1, 
an upper limit cut-off was applied to the SDSS luminosi¬ 
ties of our galaxi es, to ensure the success of the automatic 
de-blending fsee iMcConnachie et al.ll2009l f. This had as a 
result to exclude the lower-z bright-massive galaxies from 
our sample. To overcome these two biases, while comparing 
galaxies at different redshifts, we chose objects within the 
mass range of 10.4^1og(MstaT-)^11.3 Mq. This mass selec¬ 
tion can also ensure that we are not including dwarf galaxies 
in our study, which could affect our statistics. Indeed, the 
lower mass limit of log(Msteiiar) = 10.4 M© refers to SDSS 
z-band magnitudes brighter than -21.5, which correspond to 
bright galaxies, accor ding to the SDSS z- band luminosity 
functions presented in iBlanton et al.l ll200ll ~). 

Moreover, for the purposes of comparison we have also 
separated our sample into four different redshift bins, cho¬ 
sen to contain a similar number of galaxies after the mass 
selection was applied (shown in Fig. |6]in different colours). 
In Table [U we present the number and the median stellar 
mass of the galaxies in each redshift bin, before and after 
the stellar mass selection. It is evident that after the mass 
segregation, we are able to compare similar - in a statis¬ 
tical manner - samples of galaxies. In the same table we 
also present the duration of each redshift period (in Gyr), 
and the number of LTGs and ETGs, as well as dynamically 
young and old groups found in each bin. 

In §3.1, we classified the morphologies of our galaxies 
based on the fit ted radial profiles a nd the selection crite¬ 
ria described in ISimard et al.l (l201lli . Nevertheless, due to 
their larger distances, LTGs at higher-z’s might appear more 
“spheroidal” with result to be misclassified as ETGs. In 
that case, we should expect to observe more ETGs as we 
are probing at higher-z. To examine this, we have applied 
- i n addition to our c urre nt - the classification s proposed 
bv iBalogh et al.l (l2004h and iTanaka et all (l2005f) . These au¬ 
thors classified galaxies into star-forming (thus, LTG) or 
quiescent (thus, ETG), based on their Ha equivalent widths 
(EW(Ha)), with star-forming galaxies having EW(Ha)>4A. 
Comparing the results of these classifications with the one 
we performed in §3.1, we find differences of <3% at the var¬ 
ious redshift bins. In addition, using t he published galax y 
morphologies from the “Galaxy Zoo’lJ llLintott et al.l 201lll , 
we found more than 95% agreement, for classifications with 
>80% confidence, suggesting that the fraction of misclassi¬ 
fied galaxies in our sample is likely very low. 

Finally, one more probable bias associated with our 
data, could emerge from the fact that as we are probing 
at higher redshifts, the nuclear spectra will be increasingly 
contaminated by their galaxy-hosts (due to the finite size of 
the slit/fiber). Usi ng the results of the fitted radial profiles of 
our galaxies, from ISimard et ^ (1201 ill , we plot in Fig. Othe 
distributions of the ratios of the 3” SDSS spectro-fiber over 
the disk radii, at the different redshift bins. Our results show 
that at z<0.104 this contamination is insignificant, however 
at 0.104<z<0.133 around 10-15%, and at 0.133<z<0.230 
more than 20% of the disk galaxies may suffer a signifi- 


^ Available at http://www.galaxyzoo.org 


cant c ontamination in their nuclear spectra. iMoran et al.l 
1 I 2 OO 2 II showed that galaxy nuclear spectra suffering a con¬ 
tamination from their disks would move towards the lower- 
left portion of the classic BPT diagram, due to the dilution 
of their spectral lines, therefore often misclassifying AGN as 
SFN. However, a more recent study bv iMaragkoudakis et al.l 
1 I 2 OI 4 II revealed that galaxies with extra-nuclear star forma¬ 
tion can show higher [OlIl]A5007/H/3 line ratios, since lower 
metallicity Hll regions in the outer parts of galaxy discs are 
also capable of producing high-excitation emission lines, im¬ 
plying that the results of such a contamination might rather 
be challenging to interpret. 


4 RESULTS 

4.1 Nuclear activity in late-type galaxies 

As presented in the previous section, to study the evolution 
of the AGN activity in our sample, we have separated it 
into four different redshift bins. We used galaxies within the 
mass-range of 10.4^1og(Mstar)^11.3 Mq, as described in 
§3.4. The corresponding fractions of the spectroscopic clas¬ 
sifications for the different subsamples are presented in Ta¬ 
ble (2] In the discussion that follows, we will refer as AGN 
hosts, to all galaxies classified either as AGN or compos¬ 
ite objects. In Fig. [8l we present how the AGN fraction^ 
of the various sub-samples change with redshift. We notice 
that the number of AGN-hosting galaxies in our sample 
(marked with black solid line) increases by ~15% moving 
towards lower-z’s. When we separate them according to the 
dynamical state of their group (also presented in the same 
table), the AGN fractions of the LTGs found in dynami¬ 
cally old groups increase from 45% to 61%. Moreover, their 
Sy2 fractions, which is 7-8% beyond z~0.133, is changing 
to 14%-28% at lower redshifts (the range depends on the 
classification method; see §3.2). The above results suggest 
a significant increase of the AGN-hosting galaxies over the 
past 3 Gyr, accompanied by the corresponding increase of 
their Sy2-hosts fraction. On the other hand, LINER frac¬ 
tions do not seem to change at all. In contrast, the AGN 
activities of the LTGs in dynamically young groups, do not 
seem to change during the same period. Despite their AGN 
fraction seems to slightly increase (within the 1-a level), 
their Sy2 fraction remains almost constant. In Table [S] we 
present the corresponding fractions of the local isolated field 
and interacting pair LTGs samples (described in §2.2), found 
in the same mass range with our galaxies. Interestingly, we 
can notice that there are no statistical differences between 
the different environments in the local Universe, with all of 
them having AGN-hosting galaxy fractions of about 60%. 
However, we have to stress here that due to the uncertain¬ 
ties in some of the comparison samples, the results of such 
a comparison are not conclusive. 

In the first columns of Tables [2] and [3l we also present 
the number of SFN galaxies in our sample, as well as the 
comparison samples. Again, we find no differences between 
the field isolated galaxies (41%), the S-l-S and S-l-E pairs 

^ As AGN fraction, we define the fraction of AGN-hosting galax¬ 
ies in our sample (it is different from the /agn mentioned in 
§3.3) 
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Table 1. Number of galaxies and median stellar masses in the different redshift bins. 


z bin - range 

No. of galaxies“ 

\ogiMstar) [Mq]“ 

No. of galaxies^ 

log(Mstar) [Mq]'' 

T'' [Gyr] 

LTGs/ETGs'' 

DY/DO'' 

Binl 0.010-0.078 

2,237 

10.26±0.20 

1,109 

10.75±0.32 

1.0 

442/667 

61/216 

Bin2 0.078-0.104 

1,717 

10.66±0.25 

1,119 

10.80±0.32 

0.4 

359/760 

50/230 

Bins 0.104-0.133 

1,444 

10.86±0.29 

1,111 

10.85±0.33 

0.4 

330/780 

34/244 

Bin4 0.133-0.230 

1,489 

11.09±0.28 

1,016 

10.97±0.34 

1.0 

328/687 

38/217 

All redshifts 

7,417 

10.76±0.12 

4,355 

10.84±0.18 

2.8 

1,459/2,894 

183/907 


“ Total number of galaxies in each bin. 

^ Number of galaxies in each bin with stellar masses of 10.5^1og(Msta7>)^11.2 Mq. 

^ Duration of each redshift period, according to ACDM cosmology. 

^ Number of late-type (LTGs) and early-type (ETGs) galaxies in each bin, after the mass selection. 
® Number of dynamically young (DY) and old (DO) groups in each bin, after the mass selection. 


Table 2. Fractions of the type of nuclear activity in our sample, found in galaxies of 10.4^1og(Mstar)^ll-3 

Mo- 


redshift bin 

SFN“ 

Composite“ 

AGN“ 

AGN+Gomp.“ 

LINER'" 

Sy2'' 


Late-type galaxies in 

dynamically young groups (DY-LTGs) 


Binl 

39±4% 

31±4% 

30±4% 

61±3% 

13-37±4% 

16-37±5% 

Bin2 

46±6% 

33±5% 

21±3% 

54±3% 

20-22±6% 

14-18±5% 

Bin3 

43±6% 

39±6% 

19±4% 

57±3% 

25±7% 

8-18±4% 

Bin4 

48±7% 

34±6% 

18±5% 

52±3% 

15-24±6% 

12-21±6% 


Late-type galaxies in dynamically old groups (DO-LTGs) 


Binl 

39±3% 

32±3% 

29±3% 

61±3% 

21-26±4% 

14-28±4% 

Bin2 

46±4% 

35±4% 

19±3% 

54±3% 

23-29±4% 

9-25±6% 

Bin3 

50±4% 

29±3% 

21±2% 

50±2% 

16-26±4% 

24-32±6% 

Bin4 

55±5% 

30±4% 

14±3% 

45±2% 

21-29±7% 

7-8±4% 


Early-type galaxies 

in dynamically young groups (DY-ETGs) 


Binl 

23±7% 

23±7% 

54±11% 

77±6% 

41-52±12% 

15-37±7% 

Bin2 

28±7% 

25±7% 

47±10% 

72±6% 

35±11% 

15-23±7% 

Bin3 

32±11% 

37±13% 

32±11% 

68±7% 

23-31±13% 

8-23±7% 

Bin4 

45±13% 

41±13% 

14±7% 

55±5% 

8-22±8% 

0-8±8% 


Early-type galaxies in dynamically old groups (DO-ETGs) 


Binl 

28±1% 

24±1% 

48±2% 

72±2% 

33-34±3% 

13-30±2% 

Bin2 

35±2% 

26±1% 

39±2% 

65±2% 

34-38±3% 

15-18±2% 

Bin3 

39±2% 

26±1% 

35±1% 

61±2% 

30-31±3% 

13-18±2% 

Bin4 

44±2% 

26±2% 

31±2% 

56±2% 

26-30±4% 

15-22±3% 


“Classifications from the cla ssic EPT-diagram, u sing the [OIIIJA5007/H/9 vs the [NII]A6584/Ha line ratios. 
''Classifications based on the lKewlev et all 1 I 2 OO 6 II , using both the [SII](A6717-h6731)/HQ; and [OI]A6300/Ha 
lines, for galaxies already classified as AGN or Composite from the classic BPT diagram. 

Uncertainties denote the Icr Poissonian errors. 


Table 3. The nuclear activity of the comparison samples. 


Sample 

SEN 

AGN 

AGN-I-Comp. 

late-type galaxies 

Isolated field" 

41±6% 

10±3% 

59±7% 

S-l-S pairs'" 

41±7% 

34±6% 

59±8% 

S-l-E pairs'" 

42±13% 

38±13% 

58±16% 

early-type galaxies 

Isolated field" 

7±7% 

93±26% 

93±26% 

S-l-E pairs'" 

5±5% 

85±20% 

95±21% 

E+E pairs'" 

6±6% 

94±22% 

94±22% 


“Isolated field galaxies from [Hernandez-Ibaxra_e^_ah Ill2013li. with 10.4^1 og(Mstar)^11.3 M©. 

^Isolated interacting pair galaxies from iHemandez-Ibarra et al.l ll2014lh with 10.4'^log(Mstar)^ll-3 M© 
(Following these authors’ symbolism, S denotes the LTGs, and E the ETGs). 

These fractions correspond to the 87% and 82% of the field and interacting pair galaxy samples, respectively, 
which were detected to have emission lines. 

Uncertainties denote the la Poissonian errors. 
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Figure 7. Histograms of the ratios of the SDSS 3”spectro-fiber over the disk radius for the individual disk-galaxies in our sample, at 
the different redshift bins (colour-coding as in Fig. |6ll. In the upper right corners we indicate the redshift ranges and the median values 
for each bin. (A coloured version of this figure is available in the online journal) 






Figure 8. Fractions of the AGN+composite hosting galaxies over 
the total population of galaxies in each redshift bin, as a func¬ 
tion of redshift, for the four different sub-samples we presented 
in Table [ 2 ] (DY-LTGs in red, DO-LTGs in orange, DY-ETGs in 
green and DO-ETGs in blue lines), as well as for the total sam¬ 
ple (in black solid line). X-axis values are taken in the middle of 
each redshift bin, slightly displaced to resolve better. Uncertain¬ 
ties denote Itr Poissonian errors. (A coloured version of this figure 
is available in the online journal) 


(with 41% and 42%, respectively), and the dynamically 
young and old groups (both having 39% SFN at low-z’s). 
On the other hand, looking at higher-z’s, we can see that 
both dynamically young and old LTGs, have similar SFN 
(of 48% and 55%, respectively). These results are also con¬ 
sistent with our observations of their UV-optical colours. 
Studying the [NUV-r] colours of the dynamically old LTGs 
we notice significant reduction, of almost 40%, of the blue 
star forming galaxies (blue cloud) at lower redshifts, and at 
the same time an increase of 30% and 10% of the green valley 
and red sequence galaxies (quiescent galaxies), respectively 
(Bitsakis et al. - in prep.). On the contrary, this blue cloud 
decrease is not observed to be significant for the LTGs in 
dynamically young groups (only 15%). 

In a first atte mpt to interpre t these results we rely on 
the conclusions of iBitsakis et al.l (1201 ll ). for HCG galax¬ 
ies. According to these, as the time progresses the mul¬ 
tiple encounters each compact group galaxy experiences, 
will result in the built-up of its stellar mass, as well as 
to the depletion (via star formation and/or tidal strip¬ 
ping) of a significant fraction of its gas (in dynamically old 
groups is more than 70%; derived for HCG galaxies , see 
IVerdes-Montenegro et ahl I 2 OOII : iBorthakur et al.] I 2 OI 0 I) and 
dust content (more than an order of magnitude lower from 
what is o bserved in dynamica lly young and isolated field 
LTGs; see IBitsakis et aDl2014h . In dynamically old groups. 
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Table 4. Contribution of the AGN to the total IR luminosity {Jagn) of the galaxies in the 
mass range 10.4^1og(Mst(ir)^Jll-3 Mq. 


redshift bin 

/agjv 

DY-LTGs 

/agjv 

DO-LTGs 

/agjv 

DY-ETGs 

/agaj 

DO-ETGs 

Binl 

3±1% 

4±1% 

7±1% 

6±1% 

Bin2 

5±1% 

5±1% 

10±1% 

10±1% 

Bin3 

6±1% 

8±1% 

10±1% 

14±1% 

Bin4 

10±1% 

11±1% 

16±1% 

17±1% 


With DY and DO we denote the dynamically young and old groups, respectively. 
Note that for f agn the existence of an AGN is highly uncertain (see §3.3). 
Uncertainties denote the Itr Poissonian errors. 


this effect is expected to appear stronger, since they are 
denser, more dynamically evolved, and their galaxies have 
experienced most probably a higher rate of tidal encoun¬ 
ters. The high AGN fractions at lower redshifts can be ex¬ 
plained assuming that, as the star formation of the com¬ 
pact group LTGs slowly fades-out (due to the depletion of 
gas, mentioned earlier), it reveals a weaker AGN activity, 
might also been triggered by the dynamical interactions. 
This assumption is con sistent with the scenario proposed by 
iKrongold et al.l ll2002h . where galaxy interactions/merging 
initially trigger central starbursts and subsequently reveal 
a Sy2 and eventually a Syl nucleus (thus increasing the 
fraction of AGN hosting galaxies). To explain the transition 
between Sy2 and Syl nuclei, the following scenarios have 
been proposed: a very powerful, fully-obscured Syl, which 
due to feedback wil l either completely wipe-out the torus 
or make it clumpy (ISirockv et al.ll2nn8h . or a sequence on 
the AGN power that depends either on t he low accretio n 
rates or the small mass of the SMBH (e.g. lNicastroll2000l) . 
IKrongold et al.l (|2002l ) suggested that transitions from Sy2 
to Syl AGN may have been delayed as much as 1 Gyr (upper 
limit). During this period, compact group galaxies will, most 
probably, experience a large nu mber of encoun ters (having 
dynamical times of ~100 Myr; iHicksonI Il99lt) . which will 
result in the depletion of large amounts of their interstellar 
gas (via induced star formation and/or tidal stripping), thus 
preventing high accretion rates towards the SMBH, required 
to form a BLR and thus a potential S yl nuclei. This conclu¬ 
sion is consistent with the results of iNicastro et all (|2nn,3h . 
where the absence of HBLRs in Sy2 AGN can be regulated 
by the accretion rate of the SMBH. 

The above scenario predicts the existence of only a very 
small fraction of Syl nuclei in compact groups. Using the 
classihcation described in §3.1, we find that only <1%, of 
the galaxies in our sample, may have broad-lines. The almost 
complete absence of such objects does not support the incli¬ 
nation scenario, at least in most cases, according to which 
we should be a ble to observe m ore than 20% of them, in the 
local Universe iIHo et al.lll997l ). ft is possible, though, that 
heavy dust obscuration could hide them from us. Using, the 
advent of infrared extragalactic surveys, such as the WISE 
All Sky Release, we are able to unveil the mystery of the 
emission from dust obscure d regions. Bas e d on the mid-lR 
colour selection proposed bv iMateos et al.l (|2012l . see Fig. n 
we find that only 0.5% of the galaxies in our sample could 
be dust obscured AGN. In addition, according t o the WISE 
[3.4l-|' 4.6l>0.8 (Vega mags) colour selection, bv IStern et al.l 
il201^ . this fraction is also 0.3%. Therefore, it is evident 


that obscuration is not a factor that affects the absence of 
Syl sources from our sample. 

An alternative method to detect IR luminous AGN has 
b een presented in § 3.3 (and also described in more detail 
in ICiesla et al.ll2ni5h . There, we showed how our SED mod¬ 
eling is able to estimate the fraction of the total infrared 
luminosity of a galaxy, due to an AGN (/agjv). In Tabled 
we present the results of the fitting, separating galaxies ac¬ 
cording to their optical morphologies, the dynamical state 
of their group, and their redshift. Firstly, we can see the 
nearly complete absence of AGN activity in the local com¬ 
pact group galaxies, independent of morphology and group 
state (having /agn <!>%)■ On the other hand, it is also evi¬ 
dent the rise of 10% that occurs at redshifts beyond 0.133. 
We have to remin d here to the reader, that according to 
iGiesla et al.l ll2015h GIGALE’s results can be reliable only 
for fAGN>W%. Between 5%</AGiv<10%, the uncertainties 
are very high, and below 5%, its impossible to identify an 
AGN^ _ 

iMartinez et al.l (|20inh showed that even though 45% of 
HGG galaxies in their sample, host an AGN, the majority 
appear to have low-luminosities (known as low-luminosity 
AGN; LLAGN), with a median Ha luminosity of 0.71x10'*° 
erg sec“*. We estimate the corresponding value for the AGN 
hosting galaxies in our sample to be 0.76±0.01 x 10"*° erg 
sec”*, which is in excellent agreement with their result. In 
Table [S] we also present the corresponding Ha luminosities 
of the galaxies in our sample, separating them by their mor¬ 
phologies and evolutional state of their group. We can see 
that there is no statistical difference between the various 
sub-samples of LTGs and ETGs at the different redshifts, 
but there is a signihcant increase (more than an order of 
magnitude), that occurs at higher redshifts, which is ob¬ 
served in all the sub-samples. On the other hand, local in¬ 
teracting pairs and isolated galaxies display luminosities of 
~l-2xl0*° erg sec”*, which are similar to the correspond¬ 
ing compact group luminosities at higher redshifts. The fact 
that we observe a significant decrease in the AGN Ha lumi¬ 
nosities and the fAGJV, may suggest that the AGN activity 
in groups became weaker during the last 3 Gyr, possibly as 
a conseq uence of the gas starv ation scenario, initially pro¬ 
posed bv IMartinez et ^ (l2010l ). 

4.2 Nuclear activity in early-type galaxies 

In Table [2] and Fig. [S] we have also presented the corre¬ 
sponding AGN and SFN fractions of the ETGs in our com¬ 
pact group sample. Similarly to what we have seen for LTGs, 
the number of AGN-hosting ETGs found in dynamically old 
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Table 5. Median AGN Ho luminosities 


redshift bin 

DY-LTGs 

X10^®erg 

DO-LTGs 

xlO^'^erg 

DY-ETGs 

xlO^'^erg 

DO-ETGs 
x 10^*^erg 

Binl 

0.67±0.07 

0.58±0.05 

0.26±0.05 

0.34±0.02 

Bin2 

1.63±0.23 

O.SliO.lO 

0.57±0.11 

0.48±0.03 

Bin3 

1.64±0.25 

2.07±0.25 

0.76±0.21 

0.54±0.03 

Bin4 

3.70±0.63 

2.41±0.37 

1.51±0.43 

1.03±0.07 


With DY and DO we denote the dynamically young and old groups, respectively. Uncertain¬ 
ties denote the Icr Poissonian errors. 


groups is increasing towards lower redshifts, by 16%. Simi¬ 
larly, the AGN fraction of ETGs found in dynamically young 
groups seem also to increase by 15%, in the 2-a level. Ex¬ 
amining the Sy2’s in both sub-samples, it appears that in 
contrast with the dynamically young groups, where an in¬ 
crease of 10-20% is observed, in old groups these fractions 
remain high but constant. Similar is the evolution of LIN- 
ERs with the fractions of the first increasing, and those of 
the latter remaining almost constant and much higher than 
those in LTGs. Comparing our findings with those of the 
control samples (see bottom panel of Table [3} we notice 
that, emission line isolated field galaxies have a correspond¬ 
ing fraction of AGN hosts of 93%, and E-fS and E-l-E pairs 
have 95% and 94%, respectively. Unfortunately, due to the 
very large uncertainties, we cannot derive statistically sig¬ 
nificant conclusions from these comparisons. 


In an effort to interpret the observed differences be¬ 
tween the Sy2 and LINER fractions of the AGN-hosting 
ETGs in our sample, we propose the following possibilities: 
in dynamically young groups we may observe the morpho¬ 
logical transformation of Sa’s and SO’s into ellipticals. As 
discussed in §4.1, the remaining star formation activity of 
Sa’s and SO’s is expected to fade-out and subsequently reveal 
an AGN nucleus (which most probably is a Sy2). During 
that period, SO’s already at later phases of evolution may 
eventually evolv e into ellipticals (usually having LINERs; 
I Ho et all 1997bl ~l. On the other hand, by definition dynami¬ 
cally old groups host large numbers of ETGs at all epochs, 
whi ch are expected to be either evolv ed SO’s or ellipticals 
fsee lHicksorJll983 : lBitsakis et al.ll201ll L Therefore, both Sy2 
as well as LINER fractions are expected to be high at all 
times, but also constant. However, we should mention once 
more that Sy2 and LINER fractions discussed here are based 
on the AGN definition of the classic BPT diagram, while 
the true nature of LINERs is still uncertain. More specifi¬ 
cally a number of studies (e.g. lYan fc Blantonll20l3l suggest 
that the excitation mechanisms behind LINERs are a com - 
bination of post-AGB stars and shocks. ISingh et al.l (l2013l l 
argued that post-AGB stars are ubiquitous and their ion¬ 
ising effects would be present in all galaxies, having stellar 
populations older than 1 Gyr, unless they are out-shined by 
a much brighter you ng stellar population and/or a n AGN. 
On the other hand, iGonzalez-Martfn et al.l (l2006l l showed 
that at least 60% of the LINER galaxies in their sample 
could host an AGN, based on their X-ray morphologie s and 
spectra. In addition, more recently iDopita et al.l ll2015lf also 
presented the case of the LINER elliptical NGG1052, where 
an accreting SMBH was detected. 

Unfortunately, due to the morphological classification 
we applied in our sample, we cannot distinguish between 
elliptical and SO galaxies, in order to verify the proposed 
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Figure 9. AGN fractions of the compact group (in red squares 
and continuous line), isolated field (in blue triangles and dotted 
line) and pair (in orange X’s and dashed line) galaxies, as a func¬ 
tion of their stellar mass. For the purposes of comparison with the 
control samples we have selected CG galaxies at z<0.05. X-axis 
values are taken in the middle of each stellar mass bin, slightly 
displaced to resolve better. The 50% dashed line was placed to 
show where AGN rules-over SFN activity. Error-bars denote the 
la Poissonian errors. (A coloured version of this figure is available 
in the online journal) 


scenario. Nevertheless, examining their bulge-to-disk ratios 
(B/T), we can see their tendency to be more “disky” (h av- 
ing lower B/T ratios, as derived from ISimard et ^120111 ) or 
“bulgy” (with higher B/T ratios). We find that in local dy¬ 
namically old groups, ETGs have a median B/T=0.69±0.07 
that decreases to 0.62T0.06 at z>0.133. In dynamically 
young groups the corresponding B/T ratios change from 
0.64T0.18 (at Binl) to 0.52T0.25 (at z>0.133), suggesting 
that they are not significant differences in the morphologies. 

Even though the above scenarios can explain the ob¬ 
served nuclear activities for a large number of ETGs 
in our sample, there might always be the possibility 
that some of them have increased their gas content, via 
accretion/merging from gas-rich companions, which are 
more common in d ynamically young groups. Indeed, in 
iBitsakis et al.l (l2014h . we showed that about 25% of the 
Hickson compact group ETGs display blue colours and 
enhanced star formation activities. Examining their dust 
masses and star formation rates, we found that they are 
almost an order of magnitude higher in respect to those of 
field ETGs, suggesting that compact group ETGs may have 
enhanced their gas and dust content via accretion from gas 
rich members and/or the merging of dwarf companions. 
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Finally, examining the AGN Ha luminosities of the dy¬ 
namically young and old ETG samples (see Table (3), we can 
see them both decreasing more than an order of magnitude, 
as we are moving towards lower-z’s. As it was suggested for 
compact group LTGs, in §4.1, this could be an indirect in¬ 
dication of the gas depletion these sources had experienced, 
throughout time. Moreover, we can see that dynamically 
old AGN-hosting ETGs display lower Ha luminosities, at 
all redshifts, with respect to those of the dynamically young 
ETGs. These results imply that ETGs in dynamically young 
groups may have larger gas contents, either due to accretion 
from gas rich neighbours, or because they mostly host later- 
type ETGs, in the Hubble classification scheme. 


4.3 The relation between environment and AGN 
activity at lower mass galaxies 

It is generally accepted that the nuclear activity of mas¬ 
sive galaxies is mostly dominat ed by AGN r ather than cir- 
cumnucl ear star formation (i.e. Caouti[|201^^nd references 
thereinl. iHernandez-Ibarra et al] ~( 20131 7 2014h showed this 
tendency in isolated field and interacting pair galaxies. In 
figures 6 and 4 of these publications, respectively, one can 
notice that AGN in isolated field galaxies start dominating 
the emission with respect to SEN, at log(Mstar)'Sl0.7 Mq, 
wherea s in pairs this o c curs a t log(Mstar)^10.5 M©. In ad¬ 
dition, IPimbblet et al.l ll2013ll showed that AGN in clusters 
are mostly found in galaxies with log(Mstar)^10.7 Mq. In 
Fig. |9] we present the AGN fractions as a function of the 
stellar mass, for the galaxies in our sample (presented with 
the red squares and the solid line), having z< 0.05 in order to 
compare them with those of the isolated field and interact¬ 
ing pair samples (shown with blue triangles and the dotted 
line, and orange X’s and the dashed line, respectively). The 
three distributions do not appear to have any statistical dif¬ 
ference, due to their large uncertainties, yet one can see the 
tendency of the compact group galaxies to exceed the 50% 
line (where AGN rule-over SFN), at lower stellar masses, 
than the comparison samples. 

In Fig. HOI we present the same figure, this time hav¬ 
ing all compact groups found at all redshifts. We also sepa¬ 
rate CG galaxies according to the dynamical state of their 
group (into dynamically young, blue solid line, and old, red 
dashed line). It is evident that the incidence of the AGN 
activity in dynamically young groups starts prevailing over 
SFN at masses log(Mstar-)^10.6 Mq (similar to those of the 
isolated field and interacting pair galaxies), whereas in dy¬ 
namically old groups at log(Mstar)5^10.2 Mq. We perform a 
Kolmogorov-Smirnov (KS) analysis, to test the hypothesis 
the two distributions being drawn from the same parent dis¬ 
tribution. The KS probability is Pks=0.2%, thus rejecting 
the null hypothesis, and suggesting they are different, with 
a level of significance of more than 99.9%. These results sug¬ 
gest that the more dynamically evolved the group, the lower 
the stellar mass at which AGN dominate the star formation 
activity of their members. 

To interpret these results, we considered the following 
scenarios: (i) either interactions enhance AGN activity, or 
(ii) star formation fades-out faster during the evolution of 
these galaxies. Nevertheless, as we already presented there 
is no evidence to support an enhancement of the AGN ac¬ 
tivity due to interactions. On the contrary, we showed how 
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Figure 10. AGN fractions as a function of stellar mass of the 
dynamically young (in blue diamonds and continuous line) and 
dynamically old (in red triangles and dashed line) galaxies in our 
sample. Explanations as in Fig. [9] (A coloured version of this 
figure is available in the online journal) 


multiple encounters result to the decrease of the AGN lu¬ 
minosities. However, gas consnmption, du e to tidally in- 
duce d star formation (IF edotoy et a.l.l l201ll '). tidal stripping 
(i.e. [Wrdes^MoMenegro^t^iJ 20^^^_and_shock_excitation 
of th e interstellar medinm ( Cluver et al.ll201^ : lAlatalo et al.l 
l2014h . will result to the depletion of the cold gas reservoir 
of these galaxies or will make it unable to form stars (in 
the case of shocks). These results are also consistent with 
the observed absence of Syl nuclei in our sample, where the 
depletion of gas may led to small accretion rates. From the 
above it is evident that the the prevalence of AGN activity 
over star formation, is connected with the evolutional state 
of the environment (or the number of encounters galaxies 
had experienced). 


5 DISCUSSION: THE ROLE OF COMPACT 
GROUPS IN THE AGN ZOO 

The results of this study suggest that compact groups can in¬ 
fluence the evolution of the nuclear activity of their galaxies 
in an indirect manner, by affecting the available amount of 
gas. How does this result can be compared to the picture we 
already have about AGN? Do compact groups really differ¬ 
entiate their nuclear properties from galaxies in other envi¬ 
ronments? It seems that our results do not support the AGN 
unification scheme, according to which the only differentia¬ 
tion between Syl and Sy2 sources is caused by the dust ob¬ 
scuration of the torus. The almost complete absence of Syl 
nuclei in compact groups, in addition to the low AGN Ha 
luminosities, suggest that the decrease of the gas reservoir 
may le ad to low accretion r ates, which according to lNicastrd 
ll2000ll and lElitzur et al.l (l2006l l have as result to prevent 
from existence powerful Syl nuclei. The observed absence 
of such sources (also according to mid-IR colour diagnos¬ 
tics and SED modeling) is a lso consistent with what it wa s 
observed in pairs of galaxies jHernandez-dbaxra^^tQalj 201^). 
Both r esults agree with the fi ndings ofiDuItzin-Hacv an et al.l 
lll999ll and more recently iKoulouridis et ah ( 20061 1 and 
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IVillarroel fc Koriil ll20l4l , where galaxies with close compan¬ 
ions (found in distances closer than ~150 kpc; in CGs almost 
all ga laxies are found in smaller distances, see lBitsakis et all 
D) are more likely to host a Sy2 nucleus, rather than a 
Syl. 

As we m e ntion ed in §4.1 the a bove results led 
iKron gold et al.l (l2002ll . and more recentlv lKoulouridis et ahl 
(201J), to propose an evolutionary sequence in the nuclear 
activity. Initially interactions may cause the collapse of gas 
towards the inner regions of the galaxies, which will trig¬ 
ger circumnuclear star formation, as well as a weak, dust 
enshrouded AGN, which may appear as a HBLR Sy2 nu¬ 
cleus. Then, the ever increasing presence of AGN feedback 
will lead to its transformation into a Syl, and, eventually, 
the consumption and/or loss of its fuel might transform it 
back into a non-HBLR Sy2. However, the results of our cur¬ 
rent workj_as_well_as_£reviou£_AGI^om£act group studies 
(i.e. iMartinez et al.ll20inl : ISohn et al.ll201^ . suggested that 
most compact group galaxies do not seem to complete that 
process. The absence of Syl nuclei and the low AGN Ha 
luminosities, which are reducing towards higher redshifts, 
imply that this evolutionary sequence stops during its early 
phase. In §4.1, we mentioned that in typical time scales, 
where a Sy2 can turn into a Syl, galaxies in compact groups 
will most probably experience a large number of encounters, 
which will, eventually, prevent this transformation. Galaxy 
mergers also do not seem to be a possible alternative to en¬ 
rich galaxies and fuel a powerful AGN accretion. It has been 
shown that, most such events in CGs take place under dry 
cond itions, thus no gas enrichment o f the AGN-host can oc¬ 
cur JCoziol fc Plauchu-FraunI l2007l : iKonstantopoulos et al.1 

I 2 OIOII . 

Another result of this work is that the number of com¬ 
pact group AGN-hosts (especially those of galaxies found in 
dynamically old groups) increases at lower redshifts. At a 
first glance, these results seem to be inconsistent with re¬ 
cent observations, where the AGN numb er densities have 
decre ased significantly since z=l-2 (e.g. iBarger &: Cowiel 
[200i). Nevertheless, upon careful examination we find that 
it is not the number of compact group AGN-hosts that in¬ 
creases, but the fraction of SFN that decreases due to the 
depletion of gas, which affects star formation. These results 
are also confirmed by Bitsakis et al. - in prep., where the 
star formation activity of the late-type galaxies is signif¬ 
icantly decreasing since z~0.2. Moreover, in a number of 
compact group studies it is already presented how the star 
formation activity of compact group late-type galaxies can 
be seriously affected by the depletion of gas (i.e. the spe¬ 
cific star formation rates of dynamically old group LTGs 
are more than an order of ma gnitude lower to thos e of sim¬ 
ilar galaxies in the field; see iBitsakis et al.l I 2 OIII and ref¬ 
erences therein). The past interaction-triggered star forma¬ 
tion activity, the tidal stripping - r am pressure effects hav e 
been shown to be weak in HGGs (Rasmussen et al.ll2008ll . 
even though [Pesiardins et al.l (12014 1 still advocate their in¬ 
fluence - and the shock excitation of their remaining in - 
terstellar medium llGluver et al.1 [201^ : lAlatalo et al]|2014h . 
resulted either in the loss of their gas reservoir or its inabil¬ 
ity (when shocked) to form stars. These effects are expected 
to influence stronger the star formation activity that also 
occurrs in their disks, thus making them more vulnerable 
to gravitational interactions (disk material can be stripped 


easier). Moreover, according to the AGN “downsizing” sce¬ 
nario, there is an observed decrease in the AGN luminosi¬ 
ties, since z=l, which can be either explained by the de¬ 
clining accretion of the cent ral SMBH in all galax ies (the 
so-called mass-starvation; i.e. lBarger fc Cowiell2005fl . or due 
to the starvation of the massive galaxies, as the star forma- 
tion becomes dom inated by the lower stellar mass ones (i.e. 
iGowie et al.ll 199^ 1. Yet, a decrease in the AGN luminosities 
(Iagat) was also what we observed in the galaxies of our 
sample for a narrow range in stellar mass (see Table [5]), thus 
confirming that even though their AGN number densities 
increase, they seem follow the same evolution as galaxies in 
other environments. 


6 SUMMARY 

In this paper we have presented the first study of the evo¬ 
lution of the nuclear activity in compact groups of galax¬ 
ies over the past 3 Gyr (z=0.01-0.23). Our sample is also 
the largest multi-wavelength compact group sample to-date, 
consisted of 1,770 isolated compact groups (selected using 
Hickson’s criteria), containing 7,417 galaxies. We classified 
the nuclear properties of our galaxies using optical emission 
line diagrams, as well as mid-IR diagnostic methods. We 
have also fitted their UV-to-mid-IR spectral energy distribu¬ 
tions and estimated the physical properties of the galaxies, 
which were necessary for our study. Our main conclusions 
are the following: 

• We observe an increase of ~15% of the number of AGN- 
hosts found in late-type galaxies of the dynamically old 
groups, accompanied by the corresponding reduction of the 
circumnuclear star formation as we are observing galaxies 
towards lower redshifts. Examining their Sy2 galaxy frac¬ 
tions a 10-15% increase is also observed. 

• Gomparing with local isolated field and interacting pair 
galaxy samples (found at z below 0.05), we do not find dif¬ 
ferences between the various environments. Yet, this can be 
also attributed to the small number statistics of the com¬ 
parison samples. 

• Using both optical, as well as mid-IR colour classih- 
cations, we report the absence of powerful Syl nuclei (less 
than 1% for the galaxies in our sample). These results are 
also confirmed by our SED modeling. We suggest that gas 
depletion, due to past star formation activity and tidal strip¬ 
ping, can lead to low accretion rates onto the super-massive 
black hole, which in turn will result in the lack of the highly- 
accreting Syl nuclei. The decreasing accretion rates are also 
confirmed by measuring the AGN Ha luminosities of our 
galaxies, which were reduced by almost an order of magni¬ 
tude, over the past 3 Gyr. 

• An increase of about 15%, has been observed in the 
AGN fractions of the early-type galaxies in our sample. Ex¬ 
amining the observed increase of the Sy2 and LINER frac¬ 
tions in dynamically young groups, we suggest that it may 
be caused by the on-going morphological transformation of 
lenticular into elliptical galaxies. 

• We show that galaxies found in dynamically old groups, 
where they have most probably experienced a larger number 
of encounters, are more likely to host an AGN into their 
nucleus, than in dynamically young groups, at any given 
stellar mass. 
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